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Resumen 


Guajardo-Panes, R. A., Granados-Ramírez, G. R., Sánchez- 
Cohen, I, Díaz-Padilla, G., é€. Barbosa-Moreno, FE 
(septiembre-octubre, 2017). Validación espacial de datos 
climatológicos y pruebas de homogeneidad: caso Veracruz, 
México. Tecnología y Ciencias del Agua, 8(5), 157-177. 


La fiabilidad de la información climática es sustancial en 
cualquier tipo de investigación, por ello es importante 
implementar herramientas y metodologías que permitan 
identificar datos que contengan errores, ya sea por la toma 
de datos, fallas en los sensores de los instrumentos, por la 
transcripción de la información y/o por la reubicación de las 
estaciones. En el presente trabajo se describe una propuesta 
metodológica para validar e identificar registros anómalos 
de precipitación, y temperatura mínima y máxima diaria. 
La propuesta detalla el desarrollo de una metodología que 
permitió seleccionar estaciones con información mínima 
necesaria, verificar su congruencia lógica, el proceso de un 
esquema de validación espacial y la aplicación de pruebas 
de homogeneidad normal estándar (SNHT), Pettitt y 
Buishand, para validar datos climatológicos, y comprobar 
la homogeneidad de datos de precipitación y temperaturas 
registrados en estaciones del Servicio Meteorológico 
Nacional (SMN) ubicadas en el estado de Veracruz y 
áreas aledañas. Se observó que 72% de las estaciones 
con información de precipitación cumple el criterio de 
homogeneidad; en contraste, para la temperatura mínima y 
máxima, sólo el 31 y 30% de las estaciones cumplió dicho 
criterio. Lo anterior permitió fundamentar que es posible 
que los datos sean homogéneos, sin considerar ciclos o 
condiciones estacionales; sin embargo, es recomendable 


incluir el empleo de las pruebas sugeridas, tomando en 


Abstract 


Guajardo-Panes, R. A., Granados-Ramírez, G. R., Sánchez-Cohen, 
L, Díaz-Padilla, G., €: Barbosa-Moreno, F. (September-October, 
2017). Spatial validation of climatological data and homogeneity 
tests: The case of Veracruz, Mexico. Water Technology and 
Sciences (in Spanish), 8(5), 157-177. 


Given that the reliability of climate information is crucial for any 
type of research, it is important to use tools and methodologies 
that identify erroneous data resulting from data collection, faulty 
instrument sensors, transcription of information, and the relocation 
of stations. This paper describes a methodology to validate and 
identify anomalous records of rainfall and low and high daily 
temperatures. It presents: the development of a methodology to 
select stations based on the least possible amount ofinformation, and 
verify its logical coherence; the development of a spatial validation 
process; and the use of the Pettitt and Buishand Standard Normal 
Homogeneity Test (S5NHT) to validate climatological data and to 
test the homogeneity of data recorded at National Weather Service 
(SMN) stations located in the state of Veracruz and surrounding 
areas. In terms of rainfall data, we found that 72% of the stations 
met the homogeneity criterion, while for low and high temperatures 
only 31 and 30% met the criterion. Therefore, the data may be 
homogeneous, without taking into account seasonal cycles or 
conditions. Nevertheless, it is recommended that the suggested 
tests be used, taking into account cold and/or hot periods and any 
possible relationships with the dependence on rainfall (wet and/or 
dry periods). 


Keywords: Homogeneity of variance test, Standard Normal 
Homogeneity Test (5NHT), Pettitt, Buishand, rainfall, temperature 
and validation data weather. 
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cuenta los periodos de frío y /o calor, o bien la relación que 
exista con la dependencia de la precipitación (periodos 


húmedos y/o secos). 


Palabras clave: pruebas de homogeneidad de varianza, 
prueba de homogeneidad normal estándar (SNHT), Pettitt, 


Buishand, precipitación, temperatura y validación. 


Introducción 


La información climática debe ser representa- 
tiva y precisa del lugar donde las estaciones 
climatológicas se encuentran ubicadas para 
realizar estudios relacionados con la predicción 
meteorológica y climática, permitiendo que 
se reduzcan los riesgos de desastre o realizar 
predicciones a diferentes plazos para minimizar 
los riesgos de desastre que pueden generarse 
por eventos de inundación, heladas, sequías o 
el cambio climático; los registros de información 
climática deben guardar comportamientos 
homogéneos (Cao € Yan, 2012). Los datos 
climáticos pueden no representar la variación 
del tiempo actual por fallas en los instrumentos 
de medición, errores que se ocasionan de forma 
accidental por la persona responsable en la toma 
de datos, la ubicación en la que se encuentra la 
estación, entre otras, dando a lugar a variaciones 
en datos reales, y provocando que el usuario 
de información climática obtenga resultados o 
realice inferencias erróneas. Contar con series 
de datos meteorológicos de largo plazo bajo 
condiciones homogéneas resulta de interés en la 
actualidad para la comunidad científica (Costa 
á£z Soares, 2006) y por ello debe ser validada para 
desarrollar cualquier aplicación. Para alcanzar 
este objetivo es necesario aplicar metodologías 
de verificación y tratamiento que permitan iden- 
tificar las estaciones que cumplen el supuesto 
de homogeneidad en variables de temperatura 
mínima, temperatura máxima y precipitación. 


Objetivo 
En este trabajo se desarrolla una propuesta me- 


todológica para verificar que datos climáticos 
diarios registrados por estaciones del Servicio 
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Meteorológico Nacional (SMN) del estado de 
Veracruz y estados aledaños cuentan con un 
comportamiento homogéneo mediante la im- 
plementación de tres pruebas: homogeneidad 
normal estándar (SNHT), Pettitt y Buishand. 


Antecedentes 


Los investigadores con experiencia en el análisis 
de datos han citado el uso de estas pruebas que 
han empleado información climática de diversas 
latitudes del planeta, tal es el caso de Dhorde y 
Zarenistanak (2013), que utilizaron las pruebas 
de Pettitt, la prueba de homogeneidad normal 
estándar y la prueba de homogeneidad normal 
estándar desarrollada por Alexanderson y Mo- 
berg para determinar la homogeneidad en series 
de temperatura y precipitación tanto en perio- 
dos estacionales como anuales en la República 
Islámica de Irán. La aplicación de estas pruebas 
permitió clasificar las estaciones en útiles, du- 
dosas y sospechosas. Las estaciones clasificadas 
como útiles presentaron como máximo una 
inhomogeneidad en las tres pruebas aplicadas; 
las que se especificaron como dudosas fueron 
aquellas que registraron dos condiciones 
de inhomogeneidad; las catalogadas como 
sospechosas contaron con tres condiciones de 
inhomogeneidad. Se estableció como conclusión 
que las series de precipitación y temperatura 
de las 20 estaciones consideradas en el estudio 
presentaron condiciones homogéneas. Ahmad y 
Deni (2013) reportan haber empleado pruebas 
de homogeneidad normal estándar; rangos 
de Buishand, de Pettitt y de razón de Von 
Neumann para probar homogeneidad en datos 
diarios de precipitación en series registradas de 
estaciones ubicadas en la Malasia peninsular y 
detectaron que de las 59 estaciones consideradas 
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en el estudio, el 22% presentó condiciones de 
inhomogeneidad y el 78% de homogeneidad. En 
el continente europeo, Stepánek, Zahradnicek 
y Farda (2013) llevaron a cabo todo un proceso 
completo de verificación de calidad de informa- 
ción y homogeneización de registros del tiempo 
en estaciones ubicadas en la República Checa, 
aplicaron la prueba de homogeneidad normal 
estándar, la prueba bivariada de Maronna y 
Yohai, y la prueba de Easterling y Peterson. 
Estas pruebas consideraron los registros men- 
suales, estacionales y anuales de precipitación; 
el empleo de tales intervalos permitió contar 
con diversidad de resultados. En la zona del 
mar Mediterráneo, específicamente en la región 
de la Campania, al sur de Italia, Longobardi y 
Villani (2009) aplicaron las pruebas t de Stu- 
dent y una modificación de la prueba de Ward 
para comprobar la homogeneidad en series de 
precipitación en el periodo 1981-1999, a fin de 
estudiar tendencias de variación. En el sur de 
Portugal, Costa y Soares (2006) implementaron 
las pruebas de homogeneidad normal estándar, 
de rangos de Buishand y de Pettitt a datos de 
precipitación. Como resultado, 13 estaciones 
cuentan con información homogénea, mientras 
que otras 15 no cuentan con esta condición. 
En 2001, Tuomenvirta empleó la prueba de 
homogeneidad normal estándar desarrollado 
por Alexandersson (1986) para analizar tres con- 
juntos de datos: medias anuales de temperatura 
para el periodo 1961-1990 en Finlandia; series 
de temperatura y precipitación monitoreadas 
en el Atlántico Norte, y series de temperatura 
mínimas y máximas diarias correspondientes al 
periodo 1910-1995. Se encontró que los registros 
de medias anuales de temperatura y las series 
diarias de temperaturas mínimas y máximas no 
mostraron comportamientos homogéneos debi- 
do a la reubicación de estaciones, sin considerar 
el sitio original de su ubicación para minimizar 
dichos cambios. Sin embargo, para el conjunto 
de datos de precipitación y temperatura regis- 
trados en el Atlántico Norte, no se presentaron 
cambios significativos en los registros. Pérez et al. 
(2011), en Argentina, comprobaron la existencia 
de cambios abruptos en series de precipitación 


anual en el periodo 1921-2004 de 17 estaciones 
en un área de transición de clima templado hú- 
medo en el este y estepa en el oeste. Las pruebas 
empleadas fueron Buishand y Pettitt, sus resul- 
tados mostraron que se registró un aumento 
de la precipitación anual en la década de 1960 
en estaciones ubicadas en el sur de la zona de 
estudio, mientras que en el norte se presentaron 
aumentos en la década de 1970. Guentchev, 
Barsugli y Eischeid (2010) dieron a conocer una 
opción metodológica en donde se emplean las 
pruebas de Alexandersson, la prueba de rangos 
de Buishand, la prueba de Pettitt y la prueba de 
Von Neuman, para verificar tres conjuntos de 
datos de precipitación reticulares en los estudios 
de la cuenca superior del río Colorado. Los re- 
sultados indicaron que las series de tiempo son 
suficientemente homogéneas para el análisis de 
la variabilidad en el periodo 1950-1999, cuando 
se agrega en una escala subregional. Naulier 
et al. (2015) emplearon la prueba de Buishand 
para analizar las tendencias de temperatura 
en Canadá nororiental a registros dendrocro- 
nológicos y encontraron que los veranos se 
volvieron más cálidos después de 1975 y un 
incremento de los grados día durante la última 
década (2000-2010). En México, López, Sánchez 
y Vargas (2013) usaron la prueba de Pettitt, la 
de homogeneidad normal estándar (SNHT), 
Buishand y la de Von Neumann, para verificar 
la homogeneidad de los datos empleados en la 
caracterización climatológica de la cuenca baja 
del río Bravo, Tamaulipas. Por otro lado, García 
(2013) utilizó las pruebas Pettitt, homogeneidad 
normal estándar, Buishand, Von Neumamn y t 
de Student para probar homogeneidad y detec- 
tar puntos de cambio en la media para el análisis 
de la distribución de gastos máximos anuales 
en la república mexicana. Dada la importancia 
que implica contar con información climática 
confiable, el presente trabajo tiene como ob- 
jetivo describir una propuesta metodológica 
para validar datos climatológicos, y comprobar 
la homogeneidad de datos de precipitación y 
temperatura registrados en estaciones del SMN 
ubicadas en el estado de Veracruz y en áreas 
aledañas. 
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Metodología 
Delimitación geográfica del área de estudio 


Se utilizaron datos de estaciones del SMN, que 
es el organismo encargado de proporcionar 
información sobre el estado del tiempo a escala 
nacional y local de México. La información se 
obtuvo a través de la aplicación desarrollada por 
dicho organismo para Google Earth (Conagua- 
SMN, 2014), donde se despliegan en total 5 420 
estaciones climatológicas clasificadas en dos 
categorías: en operación (3 200) y suspendidas 
(2 220), las cuales se distinguen en color amarillo 
y rojo, respectivamente. La información se ex- 
trajo en formato “kmz” y se exportó a formato 
“shape” para facilitar su manejo en un Sistema 
de Información Geográfica (SIG). Se clasificaron 
las estaciones con base en su incidencia dentro 
de la cobertura vectorial de las unidades cli- 
máticas (INEGI, 2014). Dada la importancia de 
incluir en su totalidad el territorio de Veracruz, 
se delimitó un área adicional de 20 km en 
relación con el límite político del estado, para 
incluir estaciones de apoyo ubicadas en los es- 
tados vecinos de "Tamaulipas, San Luis Potosí, 
Hidalgo, Puebla, Oaxaca, Chiapas y Tabasco, así 
como considerar un área representativa de la 
precipitación en dichos estados, pues de acuer- 
do con Cruz (2013), los datos de precipitación 
son representativos en una distancia de 10 km. 


Selección de estaciones con información 
mínima necesaria 


Se emplearon registros diarios de temperatura 
máxima (TI máx), temperatura mínima (Tmín), y 
precipitación (Pp) de cada una de las estaciones; 
el periodo de tiempo considerado es desde que 
inició operaciones hasta su última fecha de 
registro. Para que la estación fuera considerada 
en el análisis, se seleccionaron aquellas que 
cumplieran los siguientes requisitos mínimos: 
1) contar con al menos 25 años de información; 
2) tener al menos 80% de registros históricos, 
y 3) que la estación no haya dejado de operar 
antes del año 2000. 
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Congruencia lógica o integridad de los datos 


Conformada la base de datos climatológica 
de las estaciones de interés, se desarrollaron 
rutinas en macros de Excel, con el compilador 
Visual Basic, para verificar los siguientes criterios 
lógicos: 1) Pp mayores o iguales a 0; 2) Tmáx 
mayores a temperaturas mínimas (Tmín); 3) 
Tmín menores a temperaturas máximas (I máx), 
y 4) verificar que los registros se registren en 
fechas congruentes. 

Los datos que no cumplieron los criterios 
fueron reemplazados con el valor -99, el cual 
se identifica como valor faltante en procesos 
posteriores. 


Validación espacial de la información 


Otro filtro empleado para verificar la informa- 
ción climática fue la validación espacial, que se 
realizó considerando lo expuesto en la norma 
UNE 500540:2004, y que describen Estévez y 
Gavilán (2008) como un nivel de validación 
opcional. El procedimiento consistió en iden- 
tificar las estaciones vecinas más cercanas a 
una estación de referencia (figura 1); a cada 
estación de referencia se calculó un intervalo de 
confianza de los registros diarios históricos de 
cada variable climática con base en la expresión: 


Há (1) 


donde x es el promedio y s la desviación 
estándar de las variables Tmáx, Tmín o Pp; 
n, el número de datos en la serie de datos cli- 
máticos, y z,, es el valor estadístico 2.58 de la 
curva normal de frecuencias asociado con un 
nivel de confianza de 99% o significancia de 1% 
(a = 0.01) (Steel € Torrie, 1985). 

Los valores registrados en las variables 
Tmáx, Tmín o Pp que no se encontraron dentro 
de los intervalos de confianza se reemplazaron 
con el valor -99 para identificarlo como valor 
faltante en procesos posteriores (figura 1 y 
cuadro 1). 

Al igual que en el apartado anterior, para 
realizar la validación espacial, se desarrollaron 
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Figura 1. Identificación de estaciones vecinas cercanas a una estación de referencia. 


rutinas en macros de Excel con el compilador 
Visual Basic; se consideraron grupos de no más 
de cinco estaciones por cada estación de referen- 
cia ubicadas dentro de un perímetro de 50 km y 
situadas dentro de la misma unidad climática. 


Aplicación de pruebas de homogeneidad 


Las pruebas empleadas en este trabajo fueron: 
1) homogeneidad normal estándar (SNHT), 2) 
Pettitt, y 3) Buishand, las cuales se describen 
brevemente a continuación, junto con sus esta- 
dísticos de prueba de acuerdo con lo planteado 
por Ahmad y Deni (2013): 


1. Homogeneidad normal estándar (SNHT), 
por sus siglas en inglés, desarrollada por 
Alexandersson (Alexandersson, 1986) y 
que modificó con Moberg (Alexandersson 
£z Moeberg, 1997); un valor T(d) compara el 
promedio de los primeros d años registrados 
con los últimos (n — d); este valor se obtiene 
con la expresión: 


T, =dZ, +(n-d)Z, para d =1,2,...n (2) 


donde: 


> (y,+7) 


1 
d S 


A Y (9:-9) (3) 
n-d S 


21 


y 2,= 


Un valor alto de T en un año d indica una 
variación “brusca”. El estadístico de prueba T, 
se define como: 


T, = máx T(d) (4) 


1sds0 


2. Pettittí prueba no paramétrica, que se 
basa en rangos y hace caso omiso de la 
normalidad de la serie. Se basa en el orden 
de rangos de los valores y. El estadístico que 
se emplea se define como: 


X, -2)1,-d(n+1) para d =1,2,...,1n (5) 


y un valor en el año m que cumple la condición: 


e 6 


es una variación fuerte en la serie (Pettitt, 1979). 


3. Buishand, por su origen bayesiano, puede 
ser empleada en variables con cualquier 
tipo de distribución y propone identificar 
un cambio en la media de la serie de datos. 
El estadístico se define como: 


d 
¡= 


S, =0 y S; => (v,-7) para d=1,2,...n (7) 
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Cuadro 1. a) Registro histórico diario de una estación de referencia de las variables Tmáx, Tmín y Pp; 


b) registros diarios de las estaciones vecinas. 


a) 
e Histórico Tmáx Tmín Pp 
Estación 
Mes Día LI LS LI LS LI LS 
30026 1 15.13 39.46 6.04 26.79 0.00 8.53 
30026 1 2 16.15 39.63 6.95 25.25 0.00 6.03 
30026 dh de le mn est 
30026 1 30 17.26 39.04 4.84 27.38 0.00 13.94 
30026 1 31 18.00 39.39 5.69 26.08 0.00 FIS 
30026 z U 17,35 39.34. 6.80 20102 0.00 0.00 
30026 2 2 16.64 40.45 6.10 25.48 0.00 0.47 
30026 ba 6do 609 609 07 bn 
30026 2 28 16.81 40.70 517 26.23 0.00 3.10 
30026 2 29 17.38 40.72 5.83 26.79 0.00 2.10 
30026 
30026 12 d 18.54 38.61 6.98 27,99 0.00 11.20 
30026 12 2 18.15 38.88 5.68 28.45 0.00 21.50 
30026 boa 669 609 609 07 2 
30026 12 30 13.58 41.47 4.67 26.97 0.00 5.05 
30026 ia 31 15.71 39.54. 4.57 26.68 0.00 216 
b) 
Estación Año Mes Día Tmáx Tmín Pp 
30075 1961 nl 1 23.50 20.00 0.00 
30075 ó0e 6ua co Eo 
30075 2009 1 1 26.00 15.50 0.00 
ES 30075 1961 1 za 20.00 19.00 0.00 
E 30075 E De sb 
5 30075 2009 1 2 27.50 15.50 0.00 
E 30075 e eN e 
= 30075 1961 il 31 23.00 16.00 0.00 
o 30075 sl se E 
o 30075 2009 1 Sl 22.00 18.00 0.00 
E 30075 det ds hos 
o 30075 1961 12 dl 26.00 18.50 0.90 
a 30075 >. sa E Ao 
E 30075 2009 12 31 27.00 16.50 0.00 
Es 30452 1956 1 1 26.50 20.00 0.00 
5 30452 0 0 eN 2 
E 30452 2007 i 1 29.00 15.00 5.50 
EE 30452 1993 1 pa 23.00 16.00 9.20 
> 30452 a o o o 
E 30452 2007 1 Za 31.50 12.50 0.00 
= 30452 1955 1 31 25.50 15.50 61.80 
Z 30452 en e A E 
E 30452 2007 1 31 30.00 19.00 0.00 
3 30452 1955 12 yl 30.00 20.00 0.00 
5 30452 de do Se 
a 30452 2005 12 31 29.00 17.00 3.50 
50 
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La prueba se define como: 


máx [Sy 
O<d=<n|s 


Q> (8) 


Otro estadístico de prueba que puede em- 
plearse es el rango que calcula la diferencia 
entre el valor mínimo y máximo de las sumas 
parciales ajustadas. La fórmula es: 


máxs; míns; 
O<sd<n O<d<n 


(9) 


Ss 


es importante agregar que Buishand propone 
valores críticos Q/vVn y R/Vn para las pruebas 
de homogeneidad (Buishand, 1982). 

El juego de hipótesis planteadas en las tres 
pruebas son las siguientes: 


H; los datos son homogéneos. 
versus 
H.; hay una fecha en la que hay un cambio en 
los datos. 


El nivel de significancia planteado es a. = 0.05, 
es decir, la probabilidad de rechazar la HA, 
cuando es verdadera (error de tipo I) es de 5% 
(Ramírez 8 López, 1993). 

Las pruebas se realizaron a las variables 
Tmáx, Tmín y Pp de cada estación con la aplica- 
ción XLSTAT en su versión de prueba, mediante 
el cual se obtuvieron los valores p (p-value) 
correspondientes; se identificaron los años en 
los que se presentó un cambio abrupto en los 
datos cuando la H, fue rechazada, y se realizó 
un concentrado de resultados para clasificar las 
estaciones en tres clases, como lo llevaron a cabo 
Schonwiese y Rapp (1997), y Wijngaard, Klein 
y Konnen (2003). 

En la clase 1 se agrupan estaciones cuyas 
pruebas de hipótesis nula fueron rechazadas a 
lo más en una de las tres pruebas, por tanto se 
considera que la información es confiable. Cuan- 
do las estaciones presentaron como máximo dos 
hipótesis nulas, las estaciones se consideraron 


en la clase 2 y la información se puede conside- 
rar como medianamente confiable. En la clase 
3 se agruparon estaciones en que se rechazaron 
tres hipótesis nulas, por lo cual la información 
se consideró como no confiable. 


Resultados 
Estaciones incluidas en el área de estudio 


En el estudio se incluyeron 260 estaciones ubi- 
cadas en el área de estudio: 162 son del estado 
de Veracruz y 98 de estados circunvecinos, las 
cuales se identifican de apoyo para procesos 
geoestadísticos y de validación. Estas estaciones 
fueron clasificadas de acuerdo con su ubicación 
en unidades climáticas, observando que 35% 
de las estaciones incide en climas subhúmedos, 
el 25% en climas húmedos y el 18% en climas 
semicálidos húmedos. Dicha clasificación se 
consideró en la validación espacial de los datos 
diarios. 


Estaciones con información mínima necesaria 


En cuanto a la verificación de información 
mínima necesaria, se observó que 226 (87%) no 
han dejado de operar antes del año 2000, con- 
tienen al menos 25 años de información y los 
registros históricos no superan el 20%. La figura 
2 muestra un panorama general de las condi- 
ciones de las estaciones climáticas en cuanto al 
inicio (línea punteada) y fin (línea continua) de 
los registros históricos (línea discontinua), y el 
total de años con información de Tmáx, Tmín y 
Pp (área sombreada). 


Verificación de congruencia lógica y 
validación espacial de la información 


Se observó que la información no congruente 
y no válida en relación con sus estaciones de 
referencia no sobrepasó el 2.83% para la varia- 
ble Tmáx; para el caso de la Tmín, no superó el 
1.54%; por último, para la variable Pp, la incon- 
gruencia y los valores no válidos en relación con 
sus estaciones vecinas no superó el 4%. 
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Figura 2. Representación del inicio y fin de registros históricos y número de años con información climática en las estaciones 


consideradas del área de estudio. 


Aplicación de pruebas de homogeneidad 


En lo correspondiente a esta etapa, se realizaron 
las pruebas de SNHT, Pettitt y Buishand a las 
variables Tmín, Tmáx y Pp de las 226 estaciones; 
con los resultados obtenidos se clasificaron las 
estaciones en confiables, medianamente confia- 
bles y no confiables. Para la variable Pp, se cla- 
sificó el 72% de las estaciones con información 
homogénea, en tanto que para la Tmín y Tmáx 
se clasificaron el 31 y 30%, respectivamente. Se 
observó que los datos correspondientes a la Pp 
presentaron condiciones de homogeneidad, en 
tanto que, contrario a lo que se esperaba, los 
datos de Tmín y Tmáx presentaron mayores 
condiciones de no homogeneidad. El cuadro 2 
muestra las estaciones desagregadas en cuanto 
a su confiabilidad por variable climática y el 
estado donde se encuentran ubicadas. 
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Con el fin de mostrar los resultados gene- 
rados de las pruebas de homogeneidad para 
las tres pruebas empleadas se conformaron tres 
tablas, donde es posible consultar los valores 
p (p-value) y los años en los que se presentó 
un cambio abrupto en las variables de Tmín 
(cuadro 4), Tmáx (cuadro 5) y Pp (cuadro 6), así 
como la clasificación de la estación con base en 
el número de H, que fueron rechazadas. 

En un sentido estricto, se seleccionaron aque- 
llas estaciones que cumplieron con la condición 
de información confiable. Se observó que para 
el estado de Veracruz, tan sólo 11 estaciones 
mostraron comportamientos homogéneos en 
las tres variables analizadas; para el estado 
de Puebla fueron seis; para San Luis Potosí se 
registraron sólo tres, y para el estado de Hidalgo 
y Oaxaca fueron 2 y 1, respectivamente (cuadro 


3 y figura 3). 
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Cuadro 2. Resultados de pruebas de homogeneidad de varianza por estado. 


Estados 
Variable | Clasificación San Luis Totales 
Chiapas | Hidalgo | Oaxaca | Puebla Pots Tabasco | Tamaulipas | Veracruz 
Confiable 4 4 16 5 1 1 40 71 
a Medianamente 
Tmín o EablE 2 1 8 1 1 20 33 
No confiable 3 11 4 15 5 2 6 76 122 
Confiable 4 3 18 6 1 1 85) 68 
A Medianamente 
Tmáx confiable 1 3 1 2 3 17 27 
No confiable 2 10 3 19 5 2 4 84 131 
Confiable 2 14 6 23 10 1 7 99 162 
Medianamente 
Pp O 1 3 5 13 22 
No confiable 3 11 1 2 1 24 42 


Cuadro 3. Estaciones con información homogénea en la zona de estudio en las variables de Tmáx, Tmín y Pp. 


Estación Estado Municipio Longitud Latitud Altitud 
13140 Jaltocán 98* 32' 18” 21” 07' 56” 232 
rr Hidalgo 

13144 Huehuetla 98 04' 34” 20* 27' 36” 466 

20294 Oaxaca San Juan Bautista Tuxtepec 96 07' 59” 18* 04” 59” 37 

21038 Guadalupe Victoria PTY A ¡OLE RO7A 2 481 

21072 Chalchicomula de Sesma 97% 27' 02” 19 01 11” 2 586 

21129 Cuyoaco 977 30" 04” IO ISOROSí 2574 

A Puebla 

21142 Venustiano Carranza 97% 40' 05” 20" 30' 37” 161 e 

21154 Esperanza 97 26' 07” SAS 2Rl 2 433 E 

21201 Hueytamalco 97% 17' 14” 19 59 14” 593 p 

24009 San Martín Chalchicuautla 98* 39 19” LO 2 187 = 

24053 San Luis Potosí | San Martín Chalchicuautla 98? 33 53” 21 30' 15” 95 2 

24113 San Vicente Tancuayalab 98% 36' 33” 217 41' 21” 40 2 

30047 Comapa 96" 41' 47” 19* 10' 20” 550 ? 

30089 Las Minas 97% 08' 51” IOSATE2SA 1 459 E 

30125 Papantla 97719 30” 20" 26' 45” 200 5 

30147 San Juan Evangelista 95 08' 45” 179 53' 00” 29) po 

30180 Chicontepec 98* 08' 28” 217 02 17” 130 A 

30327 Veracruz Las Choapas 94 09' 48” 17" 34 19” 39 > 

30345 Tantoyuca 98* 10' 39” 21" 30' 48” 94 _ 

30350 Citlaltépetl IAEZ SE 219 19 45” 210 2 

30371 Gutiérrez Zamora 97% 05' 02” 20" 26' 57" 5 S 

30377 Tamiahua 977 24' 48” AC 7 E 

30384 Jalacingo 97 18' 14” 19" 49' 42” 1749 S 
E 
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Figura 3. Distribución espacial de estaciones climatológicas consideradas en el estudio y estaciones con información 


homogénea en las variables de Tmín, Tmáx y Pp. 


Conclusiones 


Si bien la implementación de pruebas se aplicó 
a series de datos anuales, los resultados ex- 
puestos en el presente trabajo indican que los 
datos de precipitación (Pp) pueden presentar 
comportamientos homogéneos sin considerar 
ciclos o condiciones estacionales; sin embargo, 
estas condiciones deberán tomarse en cuenta 
para realizar estudios con mayor nivel de de- 
talle, es decir, para el caso en que se requieran 
estudios enfocados a una especie vegetal de 
interés, es necesario considerar ciclos de cultivo 
(primavera-verano, otoño-invierno); en caso de 
realizar comparaciones de variación del clima 
a través de distintos años, lustros y/o décadas 
para identificar variaciones climáticas deberán 
considerarse periodos estacionales tales como 


primavera, verano, otoño e invierno, o bien, 
considerar periodos de interés específico, en 
donde se requieran realizar estudios sobre 
alguno de los procesos que intervienen en el 
ciclo hidrológico, como evaporación, precipi- 
tación, transpiración, infiltración, percolación, 
afloramiento, almacenamiento y escurrimiento. 
En cuanto a los datos de temperatura mínima 
(Tmín) y temperatura máxima (TI máx), el uso 
de ciclos o condiciones estacionales puede dar 
como resultado que las series de datos muestren 
comportamientos de homogeneidad, pues se 
estarían considerando periodos con semejantes 
comportamientos de periodos de bajas y / o altas 
temperaturas. 

Aunque las pruebas que se desarrollaron 
aquí para realizar la verificación de homoge- 
neidad de las series climáticas del estado de 
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Veracruz y estados vecinos han sido empleadas 
para el mismo fin en México, una aportación 
importante del presente trabajo es la conside- 
ración de un proceso de validación espacial, 
el cual consistió en verificar registros de clima 
de cada estación con los de estaciones vecinas, 
siempre y cuando compartieran su ubicación 
con la misma unidad climática, proceso que es 
recomendable realizar previo a la aplicación de 
las pruebas de homogeneidad. Es importante 
señalar que después de realizar las pruebas de 
homogeneidad, se debe considerar un proceso 
de homogeneización de bases de datos, en 
especial a las estaciones que cuentan con datos 
medianamente confiables y no confiables, cuyo 
proceso no fue abordado en este trabajo y que 
será desarrollado en el futuro, junto con la iden- 
tificación de una prueba de las abordadas en el 
presente trabajo, a fin de identificar la de mayor 
potencia y robustez para el caso de información 
de clima en nuestro país. 

En este trabajo se muestra cómo una 
propuesta que puede ser considerada en la 
validación de datos climáticos no descarta la 
posibilidad de que sea complementaria o mejo- 
rada, incluso corregida, pues continuamente se 
cuenta con mayores conocimientos y recursos 
computacionales, que permiten el manejo de 
grandes volúmenes de datos y harán de la in- 
formación climática una base sólida y confiable 
para desarrollar trabajos de investigación con 
resultados e inferencias de mayor certidumbre. 


Agradecimientos 


Al licenciado en Sistemas Computacionales Administrati- 
vos, Moisés Fernando Cortina Cardeña, por su apoyo en el 
desarrollo de las rutinas de programación para realizar la 
validación espacial de información climática. A la licenciada 
en Estadística, Columba Falfán Castillo, por el apoyo en el 


análisis de homogeneidad a las series de datos climáticos. 


Referencias 


Ahmad, N. H,, € Deni, S. M. (2013). Homogeneity test on 
daily rainfall series for Malaysis. Matematika, 29, 141-150. 

Alexandersson, H. (1986). A homogeneity test applied to 
precipitation data. Journal of Climate, 6, 661-675. 


« ISSN 0187-8336 


Alexandersson, H., € Moeberg, A. (1997). Homogeneization 
of Swedish temperature data. Part I: Homogeneity test for 
linear trends. International Journal of Climatology, 17, 25-34. 

Buishand, T. A. (1982). Some methods for testing the 
homogeneity of rainfall records. Journal of Hydrology, 58, 
11-27. 

Cao, L. J., £ Yan, Z. W. (2012). Progress in research on 
homogenization of climate data. Advances in Climate 
Change Research, 3, 59-67. 

Conagua-SMN (2014). Estaciones climatológicas en Google 
Earth. México, DF: Comisión Nacional del Agua, 
Servicio Meteorológico Nacional. Recuperado de 
http: / /smn.cna.gob.mx /index.php?option=com_cont 
entérview=article$tid=42:normales-climatologicas-por- 
estacionércatid=16:generalé-Itemid=75. 

Costa, A. C., € Soares, A. (2006). Identification in 
inhomogeneities in precipitation time series using SUR 
models and the Ellipse test (pp. 419-428). In: Proceedings 
of Accuracy 2006. Caetano, M., €e Painho, M. (eds.). 7% 
International Syposium on Spatial Accuracy Assessment 
in Natural Resources and Environmental Sciences. 
Instituto Geográfico Portugués. Lisbon. 

Cruz, D. A. (2013). Adaptación al cambio climático en el 
área protegida Sierra Gorda, Guanajuato. México, DF: 
Universidad Nacional Autónoma de México. 

Dhorde, A. G,, € Zarenistanak, M. (2013). Three-way 
approach to test data homogeneity: An analysis of 
temperature and precipitation series over Southwestern 
Islamic Republic of Iran. Indian Geophysical Union, 17, 233- 
242. 

Estévez, J., £ Gavilán, P. (2008). Procedimientos de validación 
de datos de estaciones meteorológicas automáticas. Aplicación 
a la red de información agroclimática de Andalucía. Plataforma 
de asesoramiento y transferencia del conocimiento agrario y 
pesquero de Andalucía (pp. 1-12). Murcia: Instituto Murciano 
de Investigación y Desarrollo Agrario y Alimentario. 

García, E. (2013). Análisis de la distribución de gastos máximos 
anuales (197 pp.). Tesis de Maestría. Programa de Maestría 
y Doctorado en Ingeniería, Ingeniería Civil-Hidráulica. 
México, DF: Universidad Nacional Autónoma de México. 

Guentchev, G., Barsugli, J. J., $ Eischeid, J. (2010). 
Homogeneity of gridded precipitation datasets for the 
Colorado River basin. Journal of Applied Meteorology and 
Climatology, 49, 2404-2415. 

INEGI (2014). Información vectorial de unidades climáticas, 
escala 1:1 000 000. Aguascalientes: Instituto Nacional de 
Estadística y Geografía. 

Longobardi, A., £z Villani, P. (2009). Trend analysis of annual 
and seasonal rainfall time series in the Mediterranean 
area. International Journal of Climatology, 30, 1538-1546. 

López, Y., Sánchez, G., €: Vargas, R. (2013). Caracterización 
climatológica de la cuenca baja del rio bravo en Tamaulipas (12 
pp.). 3er Congreso Nacional de Investigación en Cambio 
Climático, Tampico, Tamaulipas, México. Recuperado 


Guajardo-Panes et al., Validación espacial de datos climatológicos y pruebas de homogeneidad: caso Veracruz, México 


de http: / /cambioclimatico-tamaulipas.org/home/ 
principal.php?page=congresos. 

Naulier, M., Savard, M. M., Bégin, C., Genanaretti, E, 
Arseneault, D., Marion, J., Nicault, A., €: Bégin, Y. (2015). 
A millennial summer temperature reconstruction for 
northeastern Canada using oxygen isotopes in subfossil 
trees. Climate of the Past, 11, 521-553. 

Pérez, S., Sierra, E., López, E., Nizzero, G., Momo, EF, €z 
Massobrio, M. (2011). Abrupt changes in rainfall in the 
Eastern area of La Pampa Province, Argentina. Theoretical 
and Applied Climatology, 103, 159-165. 

Pettitt, A. N. (1979). A nonparametric approach to the change 
point detection. Applied Statistics, 28, 126-135. 

Ramírez, M. E., £ López, O. (1993). Métodos estadísticos no 
paramétricos. Texcoco, Estado de México: Universidad 
Autónoma Chapingo. 

Schonwiese, C. D., €: Rapp, J. (1997). Climate trend atlas 
of Europe based on observation 1891-1990. International 
Journal Climatology, 18, 580-598. 

Steel, R. G., € Torrie, J. H. (1985). Bicestadística: principios y 
procedimientos. Bogotá: McGraw-Hill Interamericana. 

Stepánek, P., Zahradnicek, P., €: Farda, A. (2013). Experiences 
with data quality control and homogenization of daily 
records of various meteorological elements in the Czech 
Republic in the period 1961-2010. Idójárás - Quarterly 
Journal of the Hungarian Meteorological Service, 117, 123-141. 

Tuomenvirta, H. (2001). Homogeneity adjustments of 
temperature and precipitation series-Finnish and Nordic 
data. International Journal of Climatology, 21, 495-506. 

UNE 500540 (2004). Redes de estaciones meteorológicas 
automáticas: directrices para la validación de registros 
meteorológicos procedentes de redes de estaciones automáticas. 
Validación en tiempo real. Madrid: AENOR. 

Wijngaard, J. B., Klein, M, € Konnen, G. P. (2003). 
Homogeneity of 20th century European daily temperature 
and precipitation series. International Journal Climatology, 


23, 679-692. 


Dirección institucional de los autores 
M.G. Rafael Alberto Guajardo-Panes 


Instituto Nacional de Investigaciones Forestales, Agrícolas 
y Pecuarias 

Campo Experimental Cotaxtla 

Km. 3.5 carretera Xalapa-Veracruz, Colonia Ánimas 

91190 Xalapa, Veracruz, MÉxicO 

Tel.: +52 (55) 3871 8700, ext. 87840 
guajardo.rafaelOinifap.gob.mx 


Dra. Guadalupe Rebeca Granados-Ramírez 


Universidad Nacional Autónoma de México 
Instituto de Geografía 

Circuito Exterior s/n, Ciudad Universitaria, 
Delegación Coyoacán 

04510 Ciudad de México, MÉxIcCO 

Tel.: +52 (55) 5623 0222, ext. 45477 
rebecaQigg.unam.mx 


Dr. Ignacio Sánchez-Cohen 


Instituto Nacional de Investigaciones Forestales, 
Agrícolas y Pecuarias / Centro Nacional de Investigación 
Disciplinaria en Relación Agua, Suelo, Planta, Atmósfera 
Canal Sacramento km. 6.5, Zona Industrial 4a Etapa 
35140 Gómez Palacio, Durango, México 

Tel.: +52 (55) 3871 8700, ext. 80515 
sanchez.ignacioOinifap.gob.mx 


Dr. Gabriel Díaz-Padilla 


Instituto Nacional de Investigaciones Forestales, Agrícolas 
y Pecuarias 

Campo Experimental Cotaxtla 

Km. 3.5 Carretera Xalapa-Veracruz, Colonia Ánimas 

91190 Xalapa, Veracruz, México 

Tel.: +52 (55) 3871 8700, ext. 87841 
diaz.gabrielOinifap.gob.mx 


M.C. Finlandia Barbosa-Moreno 


Instituto Nacional de Investigaciones Forestales, Agrícolas 
y Pecuarias 

Campo Experimental Valles Centrales de Oaxaca 

Melchor Ocampo 7, Santo Domingo Barrio Bajo 

68200 Villa de Etla, Oaxaca, MÉxICO 

Tel.: +52 (55) 3871 8700, ext. 86208 
barbosa.finlandiaQinifap.gob.mx 


Tecnología y Ciencias del Agua, vol. VII, núm. 5, septiembre-octubre de 2017, pp. 157-177 


ISSN 0187-8336 + 


